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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS. 

TECHNICAL MEMORANDUM NO. 435. 

TURBULENT FLOW. * 
By L. Prandtl. 

What I am about to say on the phenomena of turbulent flow 
is still far from conclusive. It concerns, rather, the first 
steps in a new path which I hope will be followed by many 
others. 

The researches on the problem of turbulence which have 
been carried on in Gottingen for about five years, have unfor- 
tunately left the hope of a thorough understanding of turbulent 
flow very small. The photographic and kinetographic pictures 
have shown us only how hopelessly complicated this flow is. 
Figs. 1-4 are photographs of water flowing in a very long deep 
rectangular trough, as taken from above with a camera carried 
along on a car in the direction of the flow. The pictures dif- 
fer greatly according to the speed of the car, but they are all 
unpleasantly complicated. I will show motion pictures of the 
same motions at the close of my lecture. Such pictures have 
hitherto been utilized only for statistical information on the 
magnitude of the velocity variations. Otherwise, we have not 

been able to learn much from. them. That which I would call the 
*" Ueber die ausgebildete Turbulenz," a lecture delivered at 
Zurich in September, J926, before the International Congress 
for Applied Mechanics. • 
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"great problem of developed turbulence/ 1 involving a thorough 
understanding and a quantitative calculation of the processes 
through which new phenomena are constantly being developed 
and also involving the determination of that mixing force 
which is produced in each indi\ r idual instance by the contest 
of new development and damping, is not likely therefore to be 
soon solved. 

It is always possible, however, if we forego a deeper un- 
derstanding of the phenomena of turbulence, to "follow theoret- 
ically various phenomena in a logicai way controlled by exper- 
iments, especially regarding the mean motion in a given turbu- 
lent flow. Even the determination of the mean velocity, as a 
function of the place, is a technically very important task. 
The first step in this direction may be so characterized that 
the apparent viscosity forces, produced by the mixing, will be 
represented in such a form that they can be introduced into 
the hydrodynamic differential equations and thus furnish dif- t 
ferential equations for the mean motion of turbulent flows.. 

Boussinesq has already undertaken this task. He originated 
the formula, now much used, which includes an " equivalent 11 A,* 
analogous to the viscosity whereby, in addition to the molecu- 
lar tension t = \x <~ etc. , we have the apparent molecular ten- 

*I took this designation from the interesting and instructive 
article "Der Has senaus taxis ch in fre;\er Lxxft und verwandte Er- 
scheinungen," by Professor Wilhelm Schmidt, Vienna (Hamburg, 
1925). Bctissinesq and subsequent^ other hydraulic engineers 
designated this quantity by € • 



N.A.C.A.. Technical Memorandum Ho. 435 



3 



t = etc. • (1) 

3y 

in which u = the periodical mean value of the velocity compo- 
nent u. This formula has the disadvantage, however, that the 
equivalent itself depends in turn on the velocity, which must 
therefore first he found- 

After many futile attempts, I finally succeeded in obtain- 
ing an expression for the apparent friction which, although only 
a rough approximation, still quite well represents the essential 
flow characteristics of slightly viscous fluids and is. therefore 
free from the abovement ioned fault. 

In order to formulate the exchange of momentum, which pro- 
duces the apparent friction, it is customary, according to the 
method of 0. Reynolds, to divide the momentary velocity into two 
components, the periodical mea,n value and the oscillation about 
this value. 

u=u + u l V = V + V 1 w = w + w 1 

In the mean value the oscillations now give an apparent state 
of tension through the components 

p u 1 s , p u ! v 1 , etc. 

in which p = density, overlined as sign of the mean value for- 
mation. These expressions are well known, and the task is now 
to obtain for them a form in which the mean flow ("fundamental 
flow") u v w occurs. 

At this point it is essential to introduce the character- 
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is tic length for the turbulent condition, which here plays a 
similar role to that of the free wave length in the kinetic' 
gas theory. It can be considered as the diameter of the -simul- 
taneously moved fluid masses, but also as the distance trav- 
ersed by such a fluid mass before it loses its individuality 
by mixing with neighboring masses. It is easily found that 
these two distances, if the Reynolds Number* is large enough, 
differ only through a constant factor (work- of overcoming re- 
sistance in penetrating foreign fluid masses = kinetic energy 
of the mass). Following the second definition, we will desig- 
nate this length as the "mixing path 11 (or "mixing distance 11 .) 
and represent it by i. If we assume that a moving fluid mass, 
situated in a current with velocity decrease crosswise to the 
current, po.ssesses a velocity equal to the mean velocity at 
the point where it originated, and that it is. shifted by the 
mixing distance I transversely to the current, its velocity 
will then differ from the mean velocity at the new point, and 
this velocity, in its first approximation, equals I when 
the mean direction of flow relative to the X axis is chosen. 
The mean oscillation u ! can therefore be put proportional to 

I The transverse motion v f can be considered as being 

o J 

produced in the manner that two fluid masses with different 
u ! , which find themselves in one another's presence, either 

come together or move farther apart. The velocities v ! thus 

*That is, the Reynolds Number obtained from the diameter of the 
fluid mass and from the relative velocity. 
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produced can therefore he put proportional to u' . The appar- 
ent molecular tension (t 1 = p u' v T ) therefore becomes 

\3 y / 

when the proportionality factors of u' and v' and also 
the correlation factor, which would come into the formation of 
the mean product, are expressed. If the sign, which must 
change with that of ||, is taken into consideration, it is 
more correctly written 



T' = P V 



du 



ay 



ay {ci) 



Comparison with the Boussinesq form shows that agreement will 
be obtained when the "equivalent" is written 



A = p f 



ay 



(3) 



This simple formula has given very good results, in spite of its 
obvious imperfections, for the Case used as the basis of its 
dexivation, namely, a flow with velocity decrease transverse 
to the current. Above all, .in agreement with hydraulic ex- 
perience, it gives resistances proportional to the square of 
the velocity, as can be easily seen from its structure. We 
must bear in mind, however, in the results ?/hich are obtained ' 
with -this formula, that it gives only a first rough approxima- 
tion. The solutions of flow problems with this viscosity fac- 
tor show somewhat unusual characteristics, when the customary 
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viscosity factor is expressed as numerically small. Instead of 



an asymptotic transition to constant velocity, a place is found 
somewhere in the finite, where the finite sharp curve tangent i- 
ally joins a horizontal straight line. Velocity maxima always 



jy _ y i\ 3/2 m y This behavior is connected with the condition 

that for 2^ = 0 the equivalent becomes zero according to our 
3y 

formula. In reality this is not exactly what happens, since 
the equivalent does not entirely disappear, due to the disturb- 
ance in the adjacent regions. ' If, for the refinement of the 
theory, an extension of the equivalent is assumed, there is 
then obtained, at the places S~,= 0, still another equivalent 
differing from zero, and correct asymptotes and maxima with 
finite curvature are again obtained. Nevertheless, many exper- 
iments show that the curvature of the velocity curve at the 
place of the maxima is especially large, from which it is to 
be concluded that the equivalent, although not zero, is still 
considerably smaller than in the vicinity, so that the simple 
formula 3 cannot be denied a certain intrinsic justification. 

Thus far we have considered only the components of the 
molecular tension, which usually play the principal role in 
hydraulic problems. It is not difficult, however, even for 
somewhat more general cases to develop an expression compati- 
ble with the symmetry characteristics, which furnishes a full , 
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tension tensor. For example, the factor 



ay 



in formula 2 



can "be, replaced by the maximum velocity decrease transverse to 

the streamline, and the second factor by ^ by the tensor 

By 

A ti + u A . Certain experiments, to which I will recur, indi- 
cate, however, that the matter is not quite so .simple. The 
equivalent itself must rather be regarded as a tensor, so that 
in the generalization of formula 2, as soon as we pass beyond 
the planar problem, we have the product of two tensors. 

The "mixing path 11 I requires further special investiga- 
tion. It may generally be considered as a $ function of the 
place, concerning which it is first to be established, that it 
must become zero on approaching a wall, because the transverse 
motions are here prevented. The observations in smooth-walled 
tubes indicate, however, /that the viscosity still has a slight 
effect. In the field of the Blasius resistance formula I must 
be put proportional with 



in order, by using formula 2, to obtain the correct relation 
of the resistance to Reynolds Number and simultaneously the 
proportionality of the speed with the seventh root of the wall 
resistance y. * 

The relations appear to be simpler in those turbulent rno- 

tions in which no walls participate as, e.g., in the mingling 
♦Proportionality with ^ follows from v 1 ~ j for a mean 

value of v» (b = width, u = mean value). 



v e/f'f 
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of flowing fluid with the surrounding non-flowing fluid and in 
the retardation of the flow behind a moving object. In these 
cases of "free turbulence, 11 it can be assumed, at least for 
sufficiently large Reynolds Numbers that, in comparable cases, 
the phenomena are always geometrically and mechanically similar 
in a cross section transverse to the longitudinal direction. 
Consequently the mixing distances, with increasing width of 
the s tream or wake, are always proportional to the width of 
the stream, whereby the consideration at the base of formula 
2 also enables the transverse velocities to be proportional to 
the mean relative velocity u with respect to the undisturbed 
fluid. This in itself plausible assumption concerning l may 
moreover be derived, if so desired, as a result of the fact 
that the retardation of the mean motion u is produced by 
the molecular tensions given in formula 2, and that, on the 
other hand, the periodical increase in the width of a section 
of the turbulent flow, etc., is proportional to the velocity 
of the transverse motions. The results of the various ways of 
considering the subject therefore agree well with one another. 

With the addition of the law of inertia for the principal- 
motion, the ratio for the increase in the width and for the 
decrease in velocity with increasing distance from the place 
of disturbance can be' foreca-st for numerous .instances. In the 
widening of the flow, there is already developed, without the 
law of inertia (or momentum), an increase of the width prop or- 
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tional to the distance from the hole? and consequently, due to 
the law of inertia, a velocity decrease inversely proportional 
to the first power of the distance for a flow of circular 
cross section or inversely proportional to the (square) root of 
the distance in the case of a flow coming out of a slot. 

In the wake of a rod or rotational body placed transverse- 
ly to the flow, the width increases proportionally to the square 
root or cubic root of the distance, while the velocity decreases 
in inverse proportion to the square root or the 2/3 power of 
the distance- In all these rules it is, moreover , assumed that 
the velocities, or the deviations of the velocity from the un- 
disturbed flow, are small in comparison with those of the 
place of disturbance. 

The spreading of a stream with circular cross section and 
the wake of a rotational body may be taken as examples. Let b 
represent the width of the stream or wake. Let I = ab, so 
that the transverse velocity v r = I ~ ~ a u. * 

From the expression 

£b~v' (a) 
dt 

is obtained, by estimating the stream 

u 4^ ~ a u 
ax 

hence b ~ a x. 

♦This result can also be obtained through a similarity calcu- 
lation. 
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The constant momentum in all the cross sections becomes 

J ~ u 2 b 2 , whence u - .const. 

9 x 

• For the wake P = U 4^ where U is the velocity of the 
at ax 

undisturbed flow and hence 

U &° ~ a u fa) 
dx 

The momentum is J ~ p U u b 2 , which equals the resistance 
(or drag) W = c w ^ F, in which F = cross section -of body 
and c w = coefficient of drag. J = W gives 

c w U F / ' \ 

which, when introduced into formula 1, becomes 

ax 



hence b ~ y^w F x 



* / n F 

and u - Uy 

The experimental proof of the rules thus obtained, in so 
far as they relate to the widening of the stream, has already 
been given* In the wake flows, the preliminary experiments 
show certain deviations, which may possibly be connected with 
too small Reynolds Numbers. The final conclusion has not yet 
been reached. 
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If, as mentioned, we assume the mixing distance I in for- 
mula 2 to he proportional to the stream width determined by the 
above-mentioned rules and also assume that it is constant over 
the whole width at one and the sane distance, we then have suf- 
ficient data to solve the hydrodyna.mic differential equations, 
supplemented by the friction member of formula 2, in the same 
manner as customary in boundary layer calculations. (The pres- 
sure differences transverse to the direction of flow are disre- 
garded, as likewise the effects of other deformation members, 
except 3^-) On this basis -my. fellow worker, Dr. Tollmien, has 
made various calculations, which will appear shortly in " Zeit- 
schrift fur angewandte Mathematik und Mechanik. n The accompany- 
ing diagrams snow some of the results. In this connection, the 
mixing, with the neighboring air, of a broad uniform air flow, 
which comes from an opening, is treated as a special previously 
unobserved form of flow. This ca.se may be a little more closely 
investigated as an example. Here also, as in the other flow 
extension problems, I is to be put proportional to x (distance 
from opening). A dependence on y is not to be assumed, so 'that 
we can put I = ex. Here the formula is u = f{r ( ), in which 
r q = y/x. On one boundary of the field u = U; * on 'the other 
boundary u = 0 (Fig. 5). We now have' 



U 3x + 7 By p 3y d cx 



rVU 

By 



3 2 u 



By s ' 

By introducing the flow function ^ = x F(Tj) ? where 

F(tj) = f f(r\) d*H, we obtain, after a short calculation, the 
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very simple differential equation" 

F F n + 2 c 2 F" F ,n - 0 
which can be solved both through F ,f =0- (u = const.) and also 
through 

F+2c 2 F M, = 0 

The latter differential equation represents a sort of damped 'os- 
cillation, of v/hich a half-oscillation represents the function 
to be considered for the mixing region. This region is joined 
tangent ially on one side by u = U and on the other side by 
u = 0 (Fig. 6). In other cases, similar differential equations 
are obtained, most of which, however, are more difficult to solve. 

After the solution, the pressure differences transverse 
to the flow direction, which were at first disregarded, can be 
calculated from the transverse accelerations and the apparent . 
transverse tensions. In the hitherto verified examples,, these 
pressure differences were found to be less than Yjo of the dy- 
namic pressure at the maximum velocity in the given- cross sec- 
tion, so that their neglect is justified. Everywhere in the ■ 
stream there, is found a slight positive pressure. In the air 
streams of the aerodynamic laboratories, this positive pressure 
has a technical importance for precision experiments, e.g., 
for the calibration of pressure gauges for air-speed measure- 
ments. For the large Gottirfgen air stream, whose mingling with 
the surrounding air has been measured, Dr. Tollmien calculated 
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a positive pressure of about 0,006 p ^ . Perhaps the following 
figures will also be of inter est- width of mixing zone 
b = 0-255 x; mixing distance I = 0.0174 x - 0.063 b; inflow 
velocity of the air (from a state of rest) into the air stream 
to replace the air carried away by the latter, v = 0.032 U. 
The agreement of the calculated with the observed curves is 
well shown by Fig. 8, where two pressure records are given, with 
the theoretical curves represented by added dash lines. 

In addition to the above-mentioned cases, there have al- 
ready been calculated the lessening of the turbulence behind 
a screen of parallel wires, whereby the velocity deviations 
diminish in an inverse ratio to the distance, and the periodic- 
al increase in the turbulent layer, which proceeds from a sepa- 
ration layer with velocity increase. The increase' ^ occurs 
with a constant velocity proportional to the strength of the 
velocity increase. The mixing distance I = c (u x - u 2 ) t. 
(The velocity curve is here given according to formula 2 simply 
by a function of the third degree of the form Ay - By 3 .) 

A theoretical determination of the hitherto only empiric- 
ally determinable magnitude c can be obtained, I believe, in 
this as in other suitable examples, from the claim that the 
layer thus developed is no longer unstable toward small oscilla- 
tions. Of course it is not easy actually to calculate the os- 
cillations of such a layer with allowance for the apparent 
friction, even of the oscillation motion. The proof of this 
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claim is that a too small o produces too small transition 
layers or streams, which are then dynamically unstable. The in- 
stability then indicates the formation % of vortices, i.e., in- 
creased mixing, widening, etc., q.e.d. (It may also be that in- 
stability is developed for a short time in the disturbance of' 
a given wave length or oscillation period, which is then, how- 
ever, replaced' by stability. In. such cases the c value actu- 
ally involved will depend on the initial disturbances.) 

We must now discuss another question, which concerns the 
more complex problems, which cannot be so readily subjected to 
theoretical treatment. Here formula 2 or its extensions can 
be employed to determine the magnitude of the mixing distance 
in dependence on the location in the fluid, for the flow es- 
tablished by the experiment. This furnishes an especially 
clear picture of the. intensity of the mixing phenomena at each 
individual point, and we receive thereby, since it concerns 
a simple length in the mixing path, a very easily presentable, 
measure of this intensity, which can be readily transferred 
from the model to the full-scale object. Thereby remarkably 
small differences in the mean value of I are manifested in the 
previously tested cases.. In smooth troughs of uniform cross 
section, as in the case of an increasing or decreasing cross 
section and as in the violent phenomena behind an obstacle and 
in the widening of streams, mixing paths are developed, which 
constitute l/S to l/lO of the depth of the water or l/2 of 
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the width of the channel or half of the actual width of the 
stream. Figs. 9-10 give the i curves for smooth widened and 
narrowed channels- Rough channels are being tested in Got tin- 
gen, but the results have not yet been published. 

One of the most Important works of the near future will 
involve the study of the friction layers on solid bodies and 
especially the conditions for the harmful separation of the 
flow from airplane wings, the. walls of diffusers (exit cones 
of wind tunnels), etc. Here also beginnings have already been 
made. 

1 In concluding, I will speak of one more group of phenomena 
which relates to the spatial turbulence problem ( in contrast 
to the periodical or symmetric rotational problem , which alone 
has thus far been discussed)- It concerns the velocity distri- 
bution in other than cylindrical tubes. I have had especially 
careful measurements made in this connection, the results of 
which were very surprising. Instead of a distribution with 
ever more rounded 11 Isotachen" (lines of like velocity), as ob- 
tained in laminar flow, there were obtained for triangular and 
rectangular channels the lines shown in Figs. 11-12, of which 
those for the rectangular channel are still more peculiar than 
for the triangular one. An article on old observations of the 
spiral motion of water in a straight channel ("Die Wasserkraft- 
laboratorien Europas" Berlin, 1926, pp. 66-67), furnished me 
the basis for a usable explanation. The water develops "second- 
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ary motions" in all- channels of uniform but not circular cross 
section in such manner that in a corner the fluid along the 
middle of the angle flows into the corner and on both sides of 
the middle it flows out of the corner. These flows, in connec- . 
tion with the usual turbulent mixing path, enable the explana- 
tion of the observed phenomena. Momentum is ever communicated 
to the corners, thus producing the great velocities there. 
Fig. 13 shows the secondary flows for the triangular and rectang- 
ular channels. It is seen how the inward flow from the wall 
develops regions of subnormal velocity at the ends of the long 
sides and also in the middle of the short sides. For the con- 
firmation of the new views, I recently had measurements made 
with another tube, as shown in Fig. 14. The results, which 
otherwise would be very surprising, confirm the explanation in 
the best manner. The currents coming from the corners develop, 
in the middle, a double eddy, which carries the water in the 
middle line toward the level fluid and carries it away again 
at the projecting corners, thus producing hypernomal velocity 
in the middle and subnormal velocity on the sides. Such sec- 
ondary currents are generated on the surface of open troughs, 
as shown by experiments. The free surface is accordingly no 
cross section through an even flow. Moreover, the photograph- 
ically determined velocity distribution on the surface agrees 
well with that found with a Pitot tube. 
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How then are the secondary currents to be explained? In 
my opinion there is no other explanation than this. The mixing 
motion is of such a nature that, alongside the to-and-f ro motion 
in the direction of the strongest velocity gradient , there is a 
still stronger oscillation, crosswise to the former and hence 
in the direction of the 11 Isotachen. 11 If this is true, it is 
explained by a simple momentum consideration that, through this 
kind of motion, forces are generated on the convex side of the 
"Isotache, 11 whose strength is proportional to the sharpness of 
the curvature. Consideration of the mean values of (u ! )~ 
u 1 v 1 , and (v ! ) 2 give results which agree with this momentum 
consideration. 

As to why the mixing motion is of this nature, that is a 
question which belongs to the previously mentioned 11 great tur- 
bulence problem," to- which unfortunately I do not have the an- 
swer. In any event, this phenomenon plainly indicates that the 
developed turbulence is an essentially three-dimensional motion. 
This circumstance seems, however, to postpone the 'solution of 
the great turbulence problem to the very distant future, be- 
cause our present mathematical resources are quite insufficient 
for three-dimensional fluid motions. The way indicated in this 
lecture vail, however, suffice in general for the arranging of 
the experimental results, so that we can dispense with a com- 
plete explanation. 

Translation by Dwight M. Miner, 

National Advisory Committee for Aeronautics. 
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Fig, 7 

Spreading of a stream. 
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Fig. 9 Velocity distri- 
bution, 
a = -3,0° e = 1.6° 
b = -3„3° f = 0.8° 
c = -1.3° g = 3.3° 
d - 0° h = 3,9° 
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Fig. 13. 
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